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A search for new physics in the acoplanar jet topology has been performed in 2.5 fb^^ of data 
from pp colhsions at =1.96 TcV, recorded by the DO detector at the Fermilab Tevatron Colhder. 
The numbers of events with exactly two acoplanar jets and missing transverse energy are in good 
agreement with the standard model expectations. The result of this search has been used to set a 
lower mass limit of 205 GeV at the 95% C.L. on the mass of a scalar leptoquark when this particle 
decays exclusively into a quark and a neutrino. In the framework of the Little Higgs model with 
T-parity, limits have also been obtained on the T-odd quark mass as a function of the T-odd photon 
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PACS numbers: 14.80.-j, 13.85.Rm 



At hadron colliders, new colored particles predicted 
by various extensions of the standard model (SM) would 
be abundantly produced if they are light enough. The 
final state with jets and missing transverse energy (J^t) 
resulting from the decay of those particles is a promising 
channel to discover physics beyond the SM. In this Letter, 
a search for new particles in the topology consisting of 
exactly two jets and Jfx is presented using 2.5 fb^^ of 
data collected at a center-of-mass energy of 1.96 TeV with 
the DO detector during Run II of the Fermilab Tevatron 
pp Collider. The result of this search has been used to 
constrain two categories of models. 

The first category corresponds to models predicting 
the existence of leptoquarks (LQ) [l[. Those are scalar 
or vector particles carrying both a lepton and a baryon 
quantum number. They are predicted by many exten- 
sions of the SM attempting to explain the apparent sym- 
metry between quarks and leptons. To satisfy experi- 
mental constraints on flavor changing neutral current in- 
teractions, leptoquarks couple only within a single gen- 
eration. Leptoquarks decay into a charged lepton and a 
quark with a branching ratio /3, or into a neutrino and 
a quark with a branching ratio 1 — (3. Pair production 
of leptoquarks assuming /3 = therefore leads only to 
a final state consisting of two neutrinos and two quarks. 
The most stringent previous limit at 95% C.L. on the 
scalar leptoquark mass of 136 GeV [2] for /3 = was ob- 
tained by the DO collaboration with 310 pb~^ of Run II 
data. The CDF collaboration also set a lower mass limit 
of llTGeV |3] with 191 pb^^ of Run II data. Those lim- 
its, as well as the results presented in this Letter, apply 
for first- and second-generation scalar leptoquarks. For 
the third-generation, tighter limits were obtained by in- 
creasing the signal sensitivity using heavy-flavor quark 
tagging j^]. 

The second category is the Little Higgs (LH) model [1] , 
which provides an interesting scenario for physics at the 
TeV scale, predicting the existence of additional gauge 
bosons, fermions, and scalar particles with masses in 
the 100 GeV - 5 TeV range. Electroweak precision con- 
straints are satisfied by introducing a discrete symmetry 
called T-parity This symmetry is constructed such 
that all the SM states are even, while most new states 
of the LH model with T-parity (LHT) are odd. In the 
LHT model, six new Dirac T-odd quarks (T-quarks or 
Q) are the partners of the left handed T-even quarks of 
the SM. In most of the parameter space, the lightest T- 
odd particle (LTP) is the so-called "heavy photon" (Ah) 
which is stable and weakly interacting. From SM pre- 
cision measurements, it is possible to set a lower mass 
limit of ~80 GeV on the mass of Ah ■ The new parti- 
cle spectrum of the LHT model has similar properties to 
spectra of supersymmetric models. The LTP, just as the 



Lightest Supersymmetric Particle in SUSY models with 
R-parity conservation, is a dark matter candidate which 
escapes undetected. There are, however, important dif- 
ferences: the new T-odd particles have the same spin as 
their SM partner; and in the LHT model, some SM states, 
for example right-handed SM fermions or gluons, have no 
partners. In the following, the mass of the T-quarks from 
the first two generations is assumed to be degenerate, and 
pair production of those four T-quarks is considered. As 
the T-odd gauge bosons other than the Ah are relatively 
heavy, T-quarks decay into a quark and Ah in most of the 
parameter space accessible at the Tevatron. It will be as- 
sumed in the following that this branching ratio is 100%. 
Pair production of T-quarks therefore leads to a final 
state with two quarks and two LTP, giving the missing 
transverse energy signature. The only direct constraint 
from collider data on the T-quark mass is the ~100GeV 
lower limit on the mass of the supersymmetric partner 
of the first two generations quarks from LEP [8j which 
can also be applied to T-quarks. Prospective studies 
have shown that the Tevatron can be sensitive to T-quark 
masses up to ~400GeV. This sensitivity is severely re- 
duced when the mass difference between the T-quarks 
and the LTP becomes small. 

The DO detector has been described in detail in 
Ref. [lo| . Tracks are reconstructed in a silicon microstrip 
tracker and a central fiber tracker (CFT), both located 
within a 2 T superconducting solenoidal magnet. The 
liquid argon and uranium calorimeter consists of three 
cryostats. The central one covers pseudorapidities 11 1 
1 77 1 < 1.1, and the two end sections extend the coverage 
up to 1 77! « 4.2. The calorimeter is designed in projec- 
tive towers of size 0.1 x 0.1 in the (77, 0) plane, where is 
the azimuthal angle in radians. The outer muon system, 
covering \ri\ < 2, consists of tracking detectors and scin- 
tillation trigger counters in front of 1.8 T iron toroids, 
followed by two similar layers after the toroids. 

Jets were reconstructed with the iterative mid- 
point cone algorithm [12] with cone radius TZ = 
y/ (A(/))-^ + (Ay)^ — 0.5 in azimuthal angle (f> and rapid- 
ity y = ^\n[{E + Pz)/{E - p^)]. The jet energy scale 
(JES) corrections were derived from the transverse mo- 
mentum balance in photon-plus-jet events. The was 
calculated from all calorimeter cells, and corrected for 
the jet energy scale and for the transverse momenta of 
reconstructed muons. 

In events from SM processes, the presence of neutri- 
nos from W or Z decay in the final state generates large 
]^T- The main irreducible SM background in this search 
for new particles is therefore the Z{-^ i/j/)-|-jets process. 
The W{-~f lv)+]ets events also exhibit the signature, 
but their contribution can be significantly reduced by 
rejecting events with an isolated electron or muon. How- 
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ever, the charged lepton can escape detection in unin- 
strumented regions of the detector, fail identification cri- 
teria, or be a tau lepton decaying hadronically. To further 
suppress that background, events containing an isolated 
high pt track are rejected. The other SM backgrounds 
for this search are the pair production of vector bosons 
{WW, WZ, ZZ) and the production of top quarks, either 
in pairs (tt) or via the electroweak interaction. Finally, 
multijet production when one or more jets are mismea- 
sured also leads to a final state with jets and ("QCD 
background" ) . 

Events from SM processes and signal events were sim- 
ulated using Monte Carlo (MC) generators and passed 
through a full GEANxS-based [isf simulation of the de- 
tector geometry and response. They were subsequently 
processed with the same reconstruction chain as the data. 
The parton distribution functions (PDFs) used in the MC 
generators are the CTEQ6L1 ^ PDFs. A data event 
from a randomly selected beam crossing was overlaid on 
each event to simulate the additional minimum bias inter- 



actions and detector noise. The ALPGEN generator [15 1 
was used to simulate W/Z+ }ets and ti production. It 
was interfaced with PYTHIA [3| for the simulation of 
initial and final state radiation (ISR/FSR) and of jet 
hadronization. Pairs of vector bosons and electroweak 
top quark production were simulated with pythia and 
COMPHEP [17|, respectively. The next-to-leading order 
(NLO) cross sections were computed with mcfm 5.1 [3]. 
The QCD background was not simulated, since it can be 
conservatively neglected in the final stage of this analysis. 

Leptoquark pair production and decays were simulated 
with PYTHIA and the CTEQ6L1 PDFs. The LQ mass 
in the MC simulation ranged from 60 to 240 GeV. The 
NLO cross sections of this process were computed from 
a program based on [19j] with a renormalization and fac- 
torization scale (^r,f) equal to the LQ mass, and using 
the CTEQ6.1M PDF sets. 

For the LHT model, it has been shown in [9j] that T- 
quark pair production and decay to qAn is very similar 
to squark pair production and decay to qxi, where Xi 
is the lightest neutralino. Signal efficiencies were there- 
fore determined using MC events generated with pythia 
corresponding to the production and decay of these su- 
persymmetric particles. It has been checked that the spin 
differences between the T-odd particles of the LHT model 
and the supersymmetric particles do not modify the sig- 
nal efficiencies. Therefore, MC simulations of such events 
were performed to cover the Q-Ah mass plane accessible 
at the Tevatron. Concerning the signal normalization, 
the cross section of first and second generation T-quark 
pair production is equal to four times the cross section 
of heavy quark pair production, if no other new particles 
predicted by the LHT model are involved in the T-quark 
production. The NLO cross sections of this signal were 
therefore calculated using mcfm 5.1, with /ij-^f equal to 
the T-quark mass, and the CTEQ6.1M PDF sets. 



The analysis strategy follows closely the "dijet" anal- 
ysis from Ref. [23|. Events were recorded using triggers 
requiring two acoplanar jets and large I^t or where 
I/It is the vector sum of the jet transverse momenta 
[I/It ~ \^)ctsPT\)- '^^^ trigger requirements evolved 
during the Run II data taking period in order to take into 
account the increasing peak instantaneous luminosity of 
the Tevatron. At the last stage of the trigger selection, 
the requirements were typically the following: (1) I^t or 
I/It greater than 30 GeV and their separation from all jets 
greater than 25°; (2) an azimuthal angle between the two 
highest Pt jets less than 170°. Offline, events where IpT 
was higher than 40 GeV were then selected. The best 
primary vertex (PVO) was defined as the vertex with the 
smallest probability to be due to a minimum bias interac- 
tion '^I'l . The longitudinal position of PVO was required 
to be less than 60 cm from the detector center to ensure 
efficient vertex reconstruction. Good jets were defined 
as jets with a fraction of energy in the electromagnetic 
layers of the calorimeter lower than 0.95. The acopla- 
narity, i.e. the azimuthal angle between the two leading 
jets, jeti and jet2, ordered by decreasing transverse mo- 
mentum, was required to be less than 165 degrees. Then, 
the two leading jets were required to be in the central re- 
gion of the detector, with |?ydot| < 0.8, where ?7dot is the 
jet pseudorapidity calculated under the assumption that 
the jet originates from the detector center. After this 
preselection, the transverse momenta of the two leading 
jets had to be higher than 35 GeV. Finally, jets were re- 
quired to originate from the best primary vertex, based 
on their associated tracks [20"|. This was accomplished 
by requiring CPFO > 0.75, where CPFO is the fraction of 
track Pt sum associated with the jet which comes from 
PVO, CPFO= X;PT'"'''(PVO)/X;PT''''''(any PV). 

At this stage, the QCD multijet background is still 
largely dominant. To further reject those events, the 
selection criteria on IpT was increased to 75 GeV. The 
requirement that the azimuthal angle between the I^t 
and the first jet, A0(^Tjeti), exceeds 90 degrees, was 
used to remove events where a jet was mismeasured and 
generating IpT aligned to that jet. Also, the minimal 
azimuthal angle lS.(j3-adn{IpT ^ anyjet) and the maximal az- 
imuthal angle A(/)inax(^T, anyjet) between jets and I^t 
directions had to be greater than 50 degrees and lower 
than 170 degrees, respectively. 

To suppress > Zz^)-|-jets events, a veto on 

events containing an isolated electron or muon with 
Pt > 10 GeV was applied. Events with an iso- 
lated track were then rejected to further reduce that 
background. Isolated tracks were required to have 
Pt > 5 GeV, to originate from PVO with DCA{z) < 5 cm 
and DCA{r) <2cm, where DCA{z) and DCA{r) are 
the positions of the projection of the distance of clos- 
est approach between the track and PVO on the beam 
direction and in the plane transverse to the beamline, 
respectively. The number of hits in the CFT used to 
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TABLE I: Number of events observed, expected from background and signal MC simulations, and signal efficiencies for Mlq = 
200 GeV at the various stages of the analysis. The QCD multijet contribution is not included in the background contribution. 
The quoted uncertainties are the combined statistical and systematic uncertainties. 





Data 


R ft plf crni 1 i"i H 

i J clk_.J\.tii- \J lALLVA. 




StiCTTiftl efficipnpv 


PtpqpI pr'i'ioTi 


9ns 0^^ 


SO 7^12 + 'i'H'iO 

/ 0^ —1— OOOKJ 


166 ± 21 


n so? ± (W! 


1st ]eadinEr iet ttt > ^."1 OeV" 




2,S S52 + 4410 


152 ± 19 


n 276 ± n o'H4 


2nd leading jet pt > 35 GeV" 


79 985 


14 538 ± 2530 


144 ± 18 


0.262 ± 0.032 


> 75 GeV 


6,509 


5, 219 ±909 


125 ±16 


0.228 ±0.028 


A0(^T,jeti) > 90° 


6,386 


5, 148 ± 897 


124 ±15 


0.226 ±0.028 


A<;/)n.in(^T,anyjet) > 50° 


3,857 


3, 453 ± 602 


93 ±12 


0.170 ±0.021 


A<;/)n,ax(-^T,anyjet) < 170° 


2,855 


2, 568 ± 448 


81 ±10 


0.147±0.018 


Isolated electron veto 


2,347 


2, 129 ±371 


79.1 ±9.8 


0.144 ±0.018 


Isolated muon veto 


2,007 


1,880 ±328 


79.1 ±9.8 


0.144 ±0.018 


Isolated track veto 


1,472 


1,398 ±244 


73.0 ±9.1 


0.133 ±0.017 


Exactly two jets 


957 


858 ±150 


49.1 ±6.1 


0.089 ±0.011 


Final Ht cut 






optimized 




Final ]^t cut 






optimized 





"First and second jets are also required to be central (Irj^ctl < 0.8), with an electromagnetic fraction below 0.95, and to have CPFO > 0.75. 



reconstruct the track was required to be at least 8. Fi- 
nally, good quality tracks were selected by requiring the 
/dof of the track-fit reconstruction to be lower than 4. 
A hollow cone with inner and outer radii of 0.06 and 0.5 
was constructed around each track that passed those cri- 
teria. If no other track with px > 0.5 GeV and the same 
quality criteria as above was found in this hollow cone, 
the track was considered isolated. The use of a hollow, 
rather than full cone also allowed rejection of tau leptons 
decaying into three charged particles. 

Events with exactly two jets with pr > 15 GeV and 
\ridet I < 2.5 in the final state were then selected. This cri- 
terion rejects a large fraction of the remaining tt events, 
and increases the signal sensitivity at large T-quark and 
leptoquark masses once large I^t a-nd are required, 
with Ht = X]jcts^'2"' where the sum is also over all jets 
with Pt > 15 GeV and |77dot| < 2.5. Table [J summa- 
rizes the number of events observed and expected from 
MC simulations at each stage of the analysis. Figure [1] 
shows comparisons between data and MC simulations: 
the distribution of the number of jets, and the ]^t and 
Ht distributions after applying all the selection criteria 
described above. 

Finally, the two final cuts on I^t and on Ht were op- 
timized for different signals by minimizing the expected 
upper limit on the cross section in the absence of sig- 
nal. To this end and also for the final limit compu- 
tation, the CLs modified frequentist method has been 
used 2^. For the leptoquark search, two benchmarks 
were defined corresponding to low {Mlq = 140 GeV) and 
high {Mlq — 200 GeV) leptoquark masses. As summa- 
rized in Tableini the optimized values were determined to 
be Ht > 150 GeV and ^t > 75 GeV for the low mass 
selection, and Ht > 300 GeV and > 125 GeV for 
the high mass selection. In the T-quark search, five Ht- 



IpT cut combinations were used to optimally scan the 
{Q,Ah) mass plane as summarized in Table HIl In all 
cases, the contribution of the QCD multijet background 
was estimated to be small enough to be conservatively 
neglected. The number of events observed are in good 
agreement with the SM expectations. 

The uncertainty coming from the JES corrections on 
the SM backgrounds and signal predictions ranges from 
5% for low Ht and ]pT cuts to 10% for high Ht and $t 
cuts. The uncertainties due to the jet energy resolution, 
to the jet track confirmation, and to jet reconstruction 
and identification efficiencies range between 2% and 4%. 
The systematic uncertainty due to the isolated track veto 
was measured to be 3%. All these uncertainties account 
for differences between data and MC simulation, both 
for signal efficiencies and background contributions. The 
trigger was found to be fully efficient for the event sam- 
ples surviving all analysis requirements within an uncer- 
tainty of 2%. The uncertainty on the luminosity mea- 
surement is 6.1% [23|. All of these uncertainties are fully 
correlated between signal and SM backgrounds. A 15% 
systematic uncertainty was set on the Vl^/Z±jets and tt 
NLO cross sections. The uncertainty on the signal ac- 
ceptance due to the PDF choice was determined to be 
6%, using the forty-eigenvector basis of the CTEQ6.1M 



PDF set [IJ. Finally, the effects of ISR/FSR on the 



signal efficiencies were studied by varying the pythia 
parameters controlling the QCD scales and the maximal 
allowed virtualities used in the simulation of the space- 
like and time-like parton showers. The uncertainty on 
the signal efficiencies was determined to be 6%. 

The nominal NLO signal cross sections, finomj were 
computed with the CTEQ6.1M PDF and for the renor- 
malization and factorization scale /ir^' = Q, where Q was 
taken to be equal to the leptoquark or T-quark mass. 
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FIG. 1: Distribution of the number of jets (a) and distributions of I^t (b) and Ht (c) before the final optimization of the cuts 
on these two quantities; for data (points with error bars), for SM background (full histogram with shaded band corresponding 
to the total uncertainty), and for signal MC (dashed histogram). The signal drawn corresponds to pair production of scalar 
leptoquarks with Mlq = 200 GeV. 



TABLE II: For each optimized event selection, information on the signal for which it was optimized {Mlq or {Mq, M^^), and 
nominal NLO cross section), lower values of Ht and I^t selection criteria, the number of events observed, the number of events 
expected from SM backgrounds, the number of events expected from signal, and the 95% C.L. signal cross section upper limit. 
The first uncertainty is statistical and the second one is systematic. 



Mlq or {Mq,Mj^J 
(GeV) 


•^nom 

(Pb) 


{Ht,]^t) 
(GeV) 


Nobs 


-^backgrd. 




0"95 

(pb) 


Leptoquark search 


140 
200 


2.38 
0.268 


(150, 75) 
(300, 125) 


353 
12 


328±llt^? 
10.6±1.7l2:o 


229±8i23 
13.7±0.6l2:o 


1.79 
0.240 


T-quark search 


(150,100) 
(250,175) 
(300,200) 
(350,200) 
(400,150) 


59.6 
3.18 
0.868 
0.242 
0.0666 


(125,75) 
(175, 100) 
(225, 125) 
(275, 150) 
(325, 175) 


566 
147 
44 
15 
7 


513 ±Utl^ 
U0±7tll 
40 ± itl 
I3.li2.ll2:? 
4.2±1.0lJ:^ 


879 ± 167tgf 

83±12tig 
25.7 ±3.4^4;? 
16.4±1.5ttl 
10.1±0.6li;5 


17.0 
2.42 
0.780 
0.169 
0.0593 



The uncertainty due to the choice of PDF was determined 
using the full set of CTEQ6.1M eigenvectors, with the 
individual uncertainties added in quadrature. The effect 
of the renormalization and factorization scale was stud- 
ied by calculating the signal cross sections for /ir^f — Q, 
fJ-iS = Q/^ and fir,f = 2 X Q. The PDF and /ir,f effects 
were added in quadrature to compute minimum, a^^, 
and maximum, (Xmax, signal cross sections. 

For the leptoquark search, Fig. [2] shows the 95% C.L. 
observed and expected upper limits on scalar leptoquark 
production cross sections. The intersection with the 
minimal NLO cross section gives a lower mass limit of 
205 GeV for (3 = 0. The corresponding expected limit is 
207 GeV. Those limits are 214 GeV and 222 GeV, respec- 
tively, for the nominal signal cross section. 

For the T-quark search. Fig. [3] shows the 95% C.L. 
excluded regions in Q~Ah mass plane assuming that the 
branching fraction of the decay Q — > qAn is 100%. The 



largest excluded T-quarks mass, 404 GeV, is obtained for 
large mass difference between the T-quarks and the LTP. 

In summary, a search for scalar leptoquarks and for 
T-quarks produced in pp collisions at ^/s =1.96 TeV has 
been performed with a 2.5 fb""'^ data sample. This search 
was conducted in events containing exclusively two jets 
and large missing transverse energy. The results are in 
good agreement with the SM background expectations, 
and 95% C.L. limits have been set on the leptoquark 
and T-quark masses. For a single-generation scalar lep- 
toquark, a lower mass limit of 205 GeV has been obtained 
for /? = 0, improving the previous limit by 69 GeV. In the 
LHT model, limits on T-quark mass were obtained as a 
function of the Ah mass assuming 100% branching ratio 
for the decay Q qAn- T-quark masses up to 404 GeV 
are excluded when the mass difference between T-quarks 
and the LTP is large. Those are the most stringent direct 
limits to date on the T-quarks mass. 
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FIG. 2: For the leptoquark search, observed (circles) and ex- 
pected (triangles) 95% C.L. upper limits on scalar leptoquark 
production cross sections. The limits obtained with the low 
mass and high mass selections are shown separately. The 
nominal production cross sections are also shown for /? = 0, 
with shaded bands corresponding to the PDF and renormal- 
ization and factorization scale uncertainties. 
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FIG. 3: For the T-quark search, expected and observed 95% 
C.L. excluded regions in the Q-Ah mass plane. The dark 
shaded region is the observed exclusion for the minimal sig- 
nal cross section hypothesis. The light shaded band shows the 
effect on the observed exclusion coming from the theoretical 
uncertainty on the signal cross section. The full and dotted 
black lines are the observed and expected limits, respectively, 
for the nominal cross section hypothesis. The hatched re- 
gion is excluded by LEP The region below the horizontal 
dashed line {M^^ < 80 GeV) is excluded by SM precision 
measurements Q. 
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